Photo induced ionization dynamics of the nitrogen vacancy defect in diamond 
investigated by single shot charge state detection 
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The nitrogen- vacancy centre (NV) has drawn much attention for over a decade, yet detailed 
knowledge of the photophysics needs to be established. Under typical conditions, the NV can have 
two stable charge states, negative (NV - ) or neutral (NV°), with photo induced interconversion of 
these two states. Here, we present detailed studies of the ionization dynamics of single NV centres in 
bulk diamond at room temperature during illumination in dependence of the excitation wavelength 
and power. We apply a recent method which allows us to directly measure the charge state of a 
single NV centre, and observe its temporal evolution. Results of this work are the steady state NV" 
population, which was found to be always < 75% for 450 to 610 nm excitation wavelength, the 
relative absorption cross-section of NV - for 540 to 610 nm, and the energy of the NV - ground state 
of 2.6 eV below the conduction band. These results will help to further understand the photo-physics 
of the NV centre. 

PACS numbers: 76.30.Mi 



I. INTRODUCTION 

The interesting features of the negatively charged 
NV - centre are high and stable fluorescence, optical spin 
state detection, long spin coherence times, and coherent 
control of surrounding nuclear spins via hypcrfinc 
interaction [1-5]. Promising applications are the use as 
fluorescence markers [6-8], nanoscale magnetic and elec- 
tric field sensing [9-12], quantum information processing 
[13-16], and fundamental studies of spin bath physics 
[17, 18]. For all these applications, understanding and 
control of the NV centre charge dynamics is important. 

The NV consists of a substitutional nitrogen atom 
and an adjacent vacancy in the diamond carbon lattice. 
Its electronic structure is built from the three carbon 
dangling bond electrons, two electrons from the nitrogen 
lone pair orbital for the neutral charge state NV°, and 
one more electron for the negative charge state NV - . 
The zero-phonon line of NV - is at 637 nm, and phonon 
assisted sidebands of the fluorescence extend up to 
800 nm. The spectrum of NV° has a zero-phonon line at 
575 nm and also exhibits strong Stokes shifted vibronic 
bands. Both charge states can be observed in diamond. 
Our current understanding of the photochromism of 
the NV is schematically shown in Fig. 1(a) [19, 20]. 
Figure 1(b) shows real-time monitoring of the charge 
state, which we use to analyse the ionization dynamics 
[21]. 

So far, the wavelength and power dependence of the 
photo induced ionization and recombination of NV - has 
only been investigated partly [6, 19-26]. Pumping NV - 
into an unknown dark state under red illumination was 
used for high resolution microscopy in [6] . By single shot 
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NMR measurements on single NV - , it was found that 
under typically used 532 nm illumination, NV - is with 
30 % probability in this dark state, which was proposed 
to be NV° [19]. Here, we will provide further evidence 
that the unknown dark state is indeed the neutral charge 
state NV°. 

The experiments were performed using a homebuilt 
confocal microscope, and measuring the NV fluores- 
cence either by an APD or a spectrometer. We used 
a pulsed supercontinuum laser (SC400-4-PP system 
from Fianium) with 40 MHz repetition rate combined 
with an acousto optic tunable filter (AOTF) to select 
the desired wavelength. The resulting spectrum had a 
width of about 4 nm around the selected wavelength. 
Additionally, a 532 nm and a 594 nm DPSS lasers 
were used. All lasers were focused onto a diffraction 
limited spot. The used samples were ultrapure Ha type 
chemical vapour deposition (CVD) diamonds. For our 
measurements we use single shot NMR [13] and single 
shot charge state detection [21]. 

The paper is organized as follows: In Sec. II, we 
show that the above mentioned dark state of NV - 
is indeed the neutral NV°, as recently proposed. In 
Sec. Ill, we investigate the photo-induced ionization and 
recombination dynamics in dependence of excitation 
wavelength and power. With this, we can determine 
the excitation spectrum of NV - and the steady state 
population of the two charge states. In Sec. IV, the 
steady state population of NV - is measured over a 
broad spectral range of the excitation wavelength. 
Together with the results from Sec. Ill, we find that 
the population is always < 75%. Additionally, we 
determine the fluorescence saturation behaviour of NV - , 
which, combined with the population of NV - , yields 
the optimal excitation wavelength. This is found to be 
around 510 to 540 nm. Finally, in Sec. V, the ionization 



2 




NV" NV° NV" NV° NV" NV 

Valence Band 



o NV" CNV° 




Time [s 



FIG. 1. (a) Schematic illustration of the photo induced ion- 
ization and recombination of NV - , showing the ground and 
excited states of both NV - and NV°. The determination of 
the ionization energy of 2.6 eV and the recombination energy 
of 2.94 eV is shown in Fig. 7. Ionization process: 1. NV - is 
excited by one photon. 2. A second photon excites one elec- 
tron from the excited state of NV - into the conduction band, 
creating NV°. The NV will end up in the NV° ground state, 
either directly after the ionization, or possibly via the excited 
state and subsequent, fast decay. Recombination process: 3. 
NV° is excited. 4. An electron is captured from the va- 
lence band, forming NV - . (b) Timetrace of the fluorescence 
of a single NV under continuous illumination with 593 nm, 
1 fiW laser light, showing a dark state of the NV (which turns 
out to be NV°). The red line is the most likely fluorescence 
level of the NV, obtained by a hidden Markov model. Here, 
the lifetimes of the two charge states are T NV - = 56.6 ms, 
T NV o — 465 ms, and the fluorescence levels are Fl NV - = 2.2 
counts/ms, Fl NV o = 0.3 counts/ms. 



and recombination energies of NV are determined. 



II. NV" "DARK STATE" IS NV° 

By illuminating NV with light of different wave- 
lengths and measuring the fluorescence spectrum, we 
find that excitation with blue light (440 nm) populates 
the neutral charge state NV°. At these wavelengths, the 
emission of NV° is increased, whereas nearly no emission 
of NV - is observable, see Fig. 2(a). As explained in 
[19], we can use single shot NMR (see [13]) on the 
NV's nitrogen nuclear spin to probe the population 
of an electronic state of the NV. This is because the 
energy splitting of the nuclear spin depends on the 
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FIG. 2. (a) Emission spectra of a single NV by illumination 
with 440 nm (yielding the spectrum of NV°) and 514 nm 
(yielding the spectrum of NV - ). Both spectra are normal- 
ized to have the same area, (b) Measurement sequence for 
(c). First, a single shot measurement determines the state 
of the nitrogen nuclear spin. Then, a variable laser pulse is 
applied. The nuclear spin state is very robust with respect to 
excitation and ionization of the NV, and survives these pro- 
cesses. The following RF pulses flips the nuclear spin, if it is 
in resonance to a nuclear spin transition, which depends on 
the electronic spin state. Finally, the state of the nuclear spin 
is measured. Comparing the initial and final measurement 
yields the probability to flip the nuclear spin. Further details 
can be found in [19]. (c) NMR spectra for green (532 nm) il- 
lumination, for a transition (mu + 1 <-> 0) in the NV - ground 
state (left) and in the unknown dark state (right). Here, no 
extra laser pulse with 532 nm was applied, because the nu- 
clear spin state measurement includes these pulses, (d) NMR 
spectra for red (625 nm) and blue (440 nm) illumination. For 
620 - 630 nm, the laser power was 20 fiW and the pulse length 
2 ms. For 432 - 456 nm, it was 20 /xW for 10 ms. The red lines 
are Lorentzian fits. Note that for this transition, there is fast 
homogeneous dephasing, and at most 0.5 spin flip probability 
can be achieved for 100 % population of the corresponding 
electronic state [19]. 
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hyperfine interaction to the NV electron spin and the 
quadrupole interaction of the nucleus, which both 
depend on the electronic wave function being different 
for NV - and NV°. The probability to flip the nuclear 
spin at a transition frequency corresponding to a certain 
electronic state depends on the population of this state. 
The NMR measurement sequence is shown in Fig. 2(b). 
Figure 2(c) shows NMR spectra for the 14 N nuclear spin 
transition in the NV" ground state and unknown dark 
state after green (532 nm) illumination, with different 
transition frequencies. We then apply blue (440 nm) and 
red (625 nm) illumination before the NMR pulse, and 
afterwards measure the probability that the nuclear spin 
is flipped at the transition frequency corresponding to 
the dark state. The resulting NMR spectra are shown in 
Fig. 2(d). We find that both red and blue illumination 
populate this state. Because blue light populates NV° 
(see above), the unknown dark state is NV°. 



III. IONIZATION DYNAMICS 
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FIG. 3. a) Power dependence of the ionization and recombi- 
nation rates for 560 nm and 593 nm illumination wavelength. 
The solid lines are fits with quadratic power dependence, 
rate(p) = ap 2 , where p is the laser power, b) Probability 
to be in NV - in dependence of the illumination power for 
560 nm and 593 nm. 



To gain further insight into the photo induced ioniza- 
tion and recombination of NV - , we perform real-time 
monitoring of the NV charge state by measuring NV 
fluorescence [21]. We use a 650 nm long pass filter 
to mainly detect emission from NV - . Especially for 
photon energies less than the zero phonon line of NV° 
(575 nm), NV° cannot be excited, and will therefore not 
fluoresce. Therefore, NV" will give a high fluorescence 
signal, whereas NV° is dark. An example time trace of 
the NV - fluorescence with 593 nm, 1 fiW illumination 
is shown in Fig. 1(b). This real-time monitoring of 
the NV charge state is only possible if the photon 
count difference from NV" and NV° is larger than 
the photon shot noise, which is the case for low power 
illumination: The ionization rate n on and recombination 
rate rR C both depend quadratically on the excitation 
power p, i.e. non^Re p 2 . The NV" fluorescence 
/ijyy-j on the other hand, depends linearly on the 
power, fljw- °c p. Effectively, the average number 
of photons we detect from NV" before it is ionized is 
photons oc flNV- l r io-a oc |, i.e. it increases linearly 
with decreasing power. To analyse the ionization 
dynamics, we use a hidden Markov model with two 
possible fluorescence levels, which yields the most likely 
state evolution, and also the corresponding lifetimes and 
transition rates [27]. 



A. Power Dependence 

First, we performed real-time monitoring of the charge 
state (cf. Fig. 1(b)) for varying excitation power from 0.5 
- 6 fiW for 560 nm and 593 nm excitation wavelength. 
Figure 3(a) shows the ionization and recombination rates 



for these illumination parameters. The curves show the 
expected quadratic behaviour for low excitation power 
for both wavelengths. In Fig. 3(b), we calculate the pop- 
ulation of NV" in dependence of the power by detailed 
balance as 



1 NV" ~ rp ,rp ' V-V 

J NV- T J NV° 

where T NV - and T NV o are the charge state lifetimes of 
NV - and NV° obtained from the fluorescence time trace 
(here, the lifetime is the inverse of the transition rate). 
The population is nearly independent of the power, be- 
cause both rates (ionization and recombination) have a 
quadratic dependence. 

B. Wavelength Dependence 

Real-time monitoring of the charge state for varying 
excitation wavelength from 540 - 610 nm (fixed excitation 
power of 1 uW) yields the bare absorption cross-section 
of NV" and the population ratio of NV" to NV°. For 
wavelengths shorter than 540 nm fluorescence of NV~ 
and NV° cannot be resolved, because the fluorescence 
contrast as well as the charge state lifetimes decrease. 
For wavelengths longer than 610 nm, the recombination 
rate from NV° to NV" becomes too small. Figures 4(a) 
and (b) show the ionization and recombination rates 
in dependence of the illumination wavelength. For the 
recombination rate shown in Fig. 4(b), we find a local 
maximum at 575 nm, which is the zero phonon line 
of NV°. This is expected, as recent low-temperature 
studies have shown that for the recombination process, 
first NV° has to be excited, and then an electron from 
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FIG. 4. Results obtained from the ionization dynamics in 
dependence of the excitation wavelength at a fixed power of 
1 fiW. a) Ionization and b) recombination rates, c) Fluores- 
cence of NV - . Because the power is fixed, this is proportional 
to the absorption cross section of NV - . d) Probability to be 
in NV - during illumination. The local maximum at 575 nm 
is at the zero phonon line of NV° . 



the valence band is captured [20], see fig. 1(a). For 
shorter wavelengths, the recombination rate increases, 
probably due to increased absorption of NV° into the 
phononic sideband of its excited state, with a maximum 
at 550 nm. The ionization rate in Fig. 4(a) shows 
two maxima at 550 nm and 590 nm. As explained in 
the following, this results mainly from the absorption 
cross-section of NV - , i.e. due to the increased excitation 
rate of NV - at these wavelengths. In the example fluo- 
rescence time trace shown in Fig. 1 (b) , high fluorescence 
indicates emission from NV - only. This means that the 
obtained fluorescence level is directly proportional to the 
excitation efficiency (absorption cross-section) of NV - . 
In contrast to this, for high power illumination the 
real-time monitoring of the charge state is not possible, 
and observed fluorescence of the NV would be a mixture 
of the two charge states NV - and NV°. The fluorescence 
excitation efficiency of NV - in Fig. 4(c) also shows two 
maxima at 550 nm and 585-590 nm, similar to the 
ionization rate. Note that the curve does not resemble 
the mirrored fluorescence spectrum of NV - shown in 
Fig. 2(a). The double peak structure could result from 
differing vibrational modes in the ground and excited 
states, or possibly from a second electronically excited 
state involved in the excitation process [28]. 

Another important point is the population of NV - 
during illumination, because this is the desired charge 
state for most applications. This can be calculates ac- 
cording to eq. (1). The result is shown in Fig. 4(d). We 
can see that Pnv- decreases for increasing wavelength 
above 540 nm, with a local maximum at 575 nm due 



to the decreased lifetime of NV°. Populating NV - 
with > 75% probability cannot be achieved at these 
wavelengths. 



IV. NV POPULATION AND SATURATION 
BEHAVIOUR 

To obtain the population of NV - for shorter wave- 
lengths than 540 nm, we use a single shot charge state 
detection method [21]. As explained above and shown in 
Fig. 1(b), low power illumination allows direct observa- 
tion of the NV's charge state. By properly choosing the 
wavelength, power and duration of a single laser pulse, 
the number of fluorescence photons counted during this 
pulse allows single shot determination of the charge state. 
The histogram of many such measurements shows two 
Poisson distributions corresponding to the fluorescence 
of NV° and NV - , see Fig. 5 b) . By applying an arbitrary 
laser pulse before performing this charge state measure- 
ment, we can measure the probability to be in NV - after 
this laser pulse, as shown in Fig. 5(a). By fitting two 
Poisson distributions to the histogram, the ratio of the 
amplitudes yields the probability P^y- to be in NV - : 

^nv- = A , , (2) 



A 



NV- 



A 



NV° 



where ^4nv- an d ^nv° are the amplitudes of the Pois- 
sonian fit functions. Note that the first laser pulse in se- 
quence Fig. 5(a) must be long enough to reach a steady 
state population of the two charge states. 

The result is shown in Fig. 5(c). For 540 nm to 570 nm, 
we obtain similar results as above (Fig. 4). Between 
510 nm to 540 nm, the probability to be in NV - is around 
75%, whereas for < 500 nm the probability decreases. 
This means that there is no particular wavelength at 
which the NV can be initialized into NV - with > 75% 
probability. 

Another interesting point is the saturation behaviour 
of the NV under high power illumination, in dependence 
of the illumination wavelength. On the one hand, this 
yields the maximal fluorescence (saturation fluorescence) 
of the NV, which should be proportional to the probabil- 
ity to be in NV - . On the other hand, it also yields the 
saturation power, i.e. how much power is needed to get 
half of the maximum fluorescence. The desired condition 
for experiments is a high maximal fluorescence combined 
with low saturation power. For this measurement, we col- 
lected saturation curves by determining the fluorescence 
in dependence of the illumination power. The function 
used to fit the measured data and get the saturation flu- 
orescence Fs and the saturation power ig is 

where F is the fluorescence and I the illumination power 
(see Appendix A). For moderate powers we are using (< 
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FIG. 5. a) Pulse sequence for measuring the probability to 
be in NV~ for arbitrary illumination. The first laser pulse 
populates NV~ and NV° with a certain probability, which 
is measured by the second pulse. A wavelength of 594 nm 
is optimum for the measurement pulse, because absorption 
of NV - is high [cf. Fig. 4(c)], leading to high NV~ fluores- 
cence, b) Histogram of photon counts during the 594 nm 
measurement pulse, after illuminating the NV with 565 nm, 
22.5 /iW for 100 ms. The solid lines are obtained from fitting 
the sum of two Poissonian distributions. The amplitudes of 
the two Poissonians is indicated by A. The blue line indi- 
cates a threshold, which can be used to determine the charge 
state of the NV after one single measurement, c) Steady state 
probability to be in NV~ in dependence of the illumination 
wavelength, obtained from the amplitudes of the Poissonian 
fits (black dots). The data above 570 nm (orange triangles) 
is taken from Fig. 4(d). 



1 mW), this function fits the data very well. The result 
of this measurement is shown in Fig. 6. For illumination 
with 520 to 550 nm, we get high saturation fluorescence 
and low saturation power. As expected, the shape of 
the saturation fluorescence looks very similar to the 
population of NV~ at low powers shown in Fig. 5(c). 
For a three level model of NV~, one would expect the 
saturation power to be inversely proportional to the 
absorption cross-section of NV - . However, this is not 
the case here, because of the ionization processes, as 
shown in Appendix A. This is why the wavelength 
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FIG. 6. Saturation fluorescence (red dots) and saturation 
power (green triangles) of the NV, obtained from saturation 
curves. The peaks at around 513 and 530 nm are probably 
errors, and could not be confirmed by separate measurements. 
The local peak at 575 nm, however, is due to the zero phonon 
line of NV°, and can also be seen in Fig. 5(c) 



dependence of the saturation power of NV - does not 
resemble the inverse of the fluorescence of NV~ shown 
in Fig. 4(c) measured at low illumination power (which 
is proportional to the NV - absorption cross-section). 



V. IONIZATION AND RECOMBINATION 
ENERGY 

In the following, we will determine the ionization and 
recombination energy of NV - in the diamond band gap. 
For the wavelengths used in Sec. Ill (540 - 610 nm), the 
ionization and recombination rates depend quadratically 
on the illumination power, because of the underlying two- 
photon processes, see Fig. 1. However, if the photon en- 
ergy is higher than the total ionization or recombination 
energy, both processes are possible with a one photon 
process. In this case, the corresponding rate has a lin- 
ear dependence on the excitation power. Therefore, we 
measure the charge state dynamics in dependence of the 
excitation power at 435 nm to 520 nm. Because direct 
observation of the dynamics is not possible at these wave- 
lengths, we use correlated single shot charge state mea- 
surements. 

By setting a threshold in Fig. 5(b), we can determine 
the charge state of the NV non-destructively depending 
on whether the photon count number is above or below 
the threshold. This allows us to do correlated charge 
state measurements, as shown in Fig. 7(a). After the 
first charge state detection, we know the current charge 
state of the NV. Then, we apply an arbitrary laser pulse 
with certain wavelength and power, and determine the 
induced charge state dynamics by varying its duration. 
Finally, a second single shot charge state measurement is 
applied. By comparing the charge state before and after 
the arbitrary laser pulse, we get the probability that this 
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FIG. 7. a) Pulse sequence for measuring the ionization and 
recombination rates for arbitrary illumination, by correlating 
the charge state before and after the laser pulse, b) Proba- 
bility to change the charge state with the laser pulse in de- 
pendence of the pulse duration, for 470 nm and 8.4 fiW, for 
initial states NV~ and NV°. The red line is an exponen- 
tial fit. c) Ionization and recombination rates in dependence 
of illumination power p. The solid lines are parabola fits 
(f(p) — a P + bp 2 ). The error bars arise from statistical fit 
errors of the measurements shown in b) . d) Linear part of the 
ionization and recombination rates (parameter a of the above 
fit function) in dependence of the illumination wavelength. 
The red line is a fit as described in the Appendix. The error 
bars are are statistical fit errors of c). 



laser pulse changed the NV's charge state. Figure 7(b) 
shows an example measurement result for 470 nm and 
8.4 fiW, in dependence of the pulse duration. As shown 
in Appendix B, we can calculate the ionization and re- 
combination rates from this charge state dynamics for a 
specific excitation wavelength and power. The power de- 
pendence of the rates is obtained by varying the power of 
the laser pulse, and measure the corresponding dynamics 
as shown in Fig. 7(b). 

Figure 7c) shows the rates in dependence of the illumi- 
nation power for some wavelengths. We see that for de- 
creasing excitation wavelength, the curves become more 
linear. To each curve, we fit the sum of a linear and of a 



f(p) =ap + bp 2 



(4) 



where p is the laser power. The parameters a and b 
are the linear resp. quadratic part of the rates. The 
linear part a is plotted in Fig. 7(d) in dependence of the 
illumination wavelength. As explained in Appendix C, 
we create a model for the ionization and recombination 
rates to account for thermal effects and the limited 
resolution of our supercontinuum laser. This allows us 
to fit the measured data, and extract the ionization and 
recombination energy of NV - . With this model, we 
obtain an ionization energy of 2.60 eV, and a recombi- 
nation energy of 2.94 eV. 



VI. SUMMARY AND CONCLUSION 

To summarize, we have performed detailed studies of 
the photo induced ionization and recombination pro- 
cess of NV - and NV°. We found that a recently re- 
ported " dark state" of NV - is indeed the neutral charge 
state NV° (Sec. II). Real-time monitoring of the charge 
state dynamics via the NV's fluorescence was performed, 
which also yielded the bare excitation spectrum of NV~ 
(Sec. III). We also determined the population of NV - 
and the fluorescence saturation behaviour for a broad 
range of the excitation wavelength (Sec. IV). The opti- 
mal excitation wavelength for NV~ is found to be around 
510 to 540 nm. Finally, we determined the ionization en- 
ergy of NV~ to be 2.60 eV, and the recombination energy 
to be 2.94 eV (Sec. V). This was possible by observing 
the transition from a two-photon ionization and recombi- 
nation process to a one-photon process, in dependence of 
the photon energy. Note that the recombination energy 
of 2.94 eV is not the energy of the ground state of NV°: 
The ground state of NV° can also be seen as a hole in 
the ground state of NV - , and the photon has to provide 
the energy difference of the valence band and the ground 
state of NV - for the recombination. This means that 
the ground state of NV - lies 2.6 eV below the conduc- 
tion band and 2.94 eV above the valence band. The sum 
of these two is 5.54 eV, close to the diamond band gap 
of 5.48 eV [29]. 

These results provide further insight into the proper- 
ties and photo-physics of the NV. A possible application 
of this process is stochastic high resolution microscopy, 
similar to PALM and STORM, which requires the acti- 
vation of only few emitters in a diffraction limited area at 
the same time. This is in principle possible with the NV, 
even in a reversible fashion, as can be seen in Fig. 1(b). 
We expect that these results are not sensitive to the envi- 
ronment of the NV, but only depend on the band struc- 
ture of the diamond. Only if the inherent ionization / 
recombination rates (in dark) are comparable to the mea- 
surement time, e.g. through shifts of the Fermi level (cf. 
[30]), an effect is expected. A remaining question is the 
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reason for the shape of the absorption spectrum of NV - 
shown in Fig. 4(c), and whether this two-maxima struc- 
ture can be explained by the transition probabilities and 
density of states according to the vibrational modes in 
the ground and excited states, or if it arises from some 
other wavelength dependent mechanism. 
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Appendix A: Saturation Behaviour 

Here, we discuss the expected saturation behaviour of 
NV~. We use a four level model, which includes the 
NV~ ground state (G), excited state (E) and metastable 
state (M), and one effective state for NV° (0), see Fig. 8. 
We do not model the dynamics within NV°, because we 
are only interested in the dynamics of NV~. The rate 
equations for these four states are 

Pg = -I&PG + AegPe + AmgPm + r le p , 

pE = Icrp G + (-A E G - A E M - Mon)PE, 

Pm = AemPe — AmgPm, 

Po = IvionPE - r rc p , (Al) 

where pi denotes the population of state i, Xij a decay 
rate from state i to state j, I the illumination intensity, a 
the absorption cross-section of NV~ , and a lon the absorp- 
tion cross-section of NV~. The rate r re is the effective 
recombination rate, which is 



(A2) 



where a vc is an effective recombination cross-section, and 
I is the saturation power of NV°. This formula is has 
two components, the population of the NV° excited state, 
similar to eq. (3), multiplied by a term proportional to 
the illumination intensity /. By setting all time deriva- 
tives in (Al) to zero, 



and using 



=p G =p E = Pm =Po, 



1 = Pg + Pe + Pm +Pa, 



(A3) 



(A4) 



we can calculate the steady state populations of the four 
states. We want to determine the expected fluorescence 




FIG. 8. Illustration of the energy levels and rates, which are 
used to model the saturation behaviour. 



of the NV, in dependence of the illumination intensity. 
The fluorescence of NV~ is proportional to the popula- 
tion of the NV - excited state E, for illumination with 
moderate laser intensity, i.e. as long as the fluorescence 
decay rate Aeg is much larger as the ionization rate /<7i on , 
which is the case for our measurements. This population 
can be written as 



Pe =Ps t , T , 
where the saturation population ps is 



PS 



1 

1 + Aem/AmG + O'ion/cTrmre + &ion/&' 



and the saturation power Is is 

Aeg + Aem + Iq^cu 



Is = 



a + CtAem/AmG + Cion + OTTion/oV 



(A5) 



(A6) 



(A7) 



Equation (A5) is proportional to the fit function eq. 3, 
justifying its use. Additionally, eq. (A7) shows that 
the saturation power is not inversely proportional to 
the absorption cross-section of NV - , as would be the 
case without the ionization process. Instead, it also 
depends on the ionization dynamics, and especially the 
saturation power I of NV° . This is probably the reason 
for the steep increase of the saturation power in Fig. 6 
at above 570 nm. 



Appendix B: Ionization and recombination rate 
equations 

The charge state dynamics of the NV can be described 
by a system of two first order linear rate equations 



Pnv° 



— A-o Ao- 
A-o — Ao- 



Pnv- 
Pnv° 



(Bl) 



p 



where Pnv-i Pnv° are the population probabilities of 
NV - resp. NV°, and A_o, Aq- are the rates from NV - 
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to NV° resp. the other direction. The time evolution of 
p is then given by 



Pit) 



At 



P(0), 



(B2) 



which can be solved by diagonalizing A. In our case, the 
starting point p(0) after the first charge state detection 

measurement is always cither ^ q (f° r NV~) or ^ ^ 

(for NV°). Then, the evolution is given by 



P-/o(t) 



A-o+An- 

A-o 
A-o+Ao- 



±- 



A °-/-° e -(A-o+A -)t ( -1 

"A_o + Ao- ' V 1 

(B3) 

where p_/ (t) indicates the evolution after initializing in 
NV~ resp. NV°. These two equations describe our mea- 
surement for the flip probabilities as shown in Fig. 7(b). 
From the steady state population p(oo) and the pumping 
rate A = A_o + Ao- provided by fitting the measurement 
results yield the two rates A_o and Ao-, which are shown 
in Fig. 7(c). 



Appendix C: Model for measured ionization and 
recombination rates 



The basic idea is that due to thermal effects and due to 
the spectral width of our supercontinuum laser (of around 
4 nm), we get a probability distribution for the energy 
absorbed by the electron (for NV~) or by the hole (for 
NV°). Since we do not know the energy distribution of 
the NV~ ground state due to thermal effects, we simply 
assume a Gaussian distribution P(Ef) for the final energy 



Ef of the electron (hole) after absorbing a photon, to 
approximately account for both the thermal effects and 
the spectral width of the photon, i.e. 



P(E t ) oc e 



(CI) 



where huj is the desired photon energy, and a is the width 
of the distribution, which will also be a fit parameter. 
The total ionization rate is given by the convolution of 
the final energy of the electron (hole) with the density 
of states (DOS) for electrons in the conduction band (for 
holes in the valence band). We approximate the band 
structure at the band edge with a parabola, such that 
the DOS g is the same as for a free electron gas, which is 



g(E) cx y/E - Eq, E > E 



(C2) 



where E is the energy and E the band edge energy rel- 
ative to the ground state of NV~ . E is the parameter 
we would like to determine. The rate therefore is 



r(hu) = A y/Ef - E e 



-2\— ) 



dE f . (C3) 



This can be solved numerically and fit to the measure- 
ment data. For the ionization rate, not all the data was 
used, because the approximation for the DOS in the band 
is only valid near the band edge. Therefore, only the data 
up to 445 nm was used. For the recombination rate, only 
four measurements points, which show a one photon pro- 
cess, are given. Therefore, the dependency of the two fit 
parameters E and a is very high. For this fit, we used 
the same value for a = 0.069 eV as obtained by the fit 
for the ionization rate. This is also physically motivated, 
because the mechanism is the same, the spectral width 
of the photon and thermal effects. 
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